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ABSTRACT 

We present a self-consistent mean field theory of the dynamo in 3D and tur- 
bulent diffusion in 2D in weakly ionized gas. We find that in 3D, the backreac- 
tion does not alter the beta effect while it suppresses the alpha effect when the 
strength of a mean magnetic field exceeds the critical value B c ~ \J v in T n (v 2 ) / R m . 
Here, z/j n , r n , and R m are ion- neutral collision frequency, correlation time of neu- 
trals, and magnetic Reynolds number, respectively. These results suggest that a 
mean field dynamo operates much more efficiently in weakly ionized gas where 
VinTn ^> 1, compared to the fully ionized gas. Furthermore, we show that in 2D, 
the turbulent diffusion is suppressed by back reaction when a mean magnetic 
field reaches the same critical strength B c , with the upper bound on turbulent 
diffusion given by its kinematic value. Astrophysical implications are discussed. 

Subject headings: ISM:magnetic fields — MHD — turbulence 

1. Introduction 

One of the most outstanding problems in the astrophysical MHD is to explain the origin 
of ubiquitous magnetic fields in stars, galaxies, interstellar medium (ISM), etc. These mag- 
netic fields are often observed to be coherent on scales much larger than the characteristic 
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scale of turbulence, with their energy being comparable with fluid kinetic energy (i.e., in 
equipartition). For instance, galactic magnetic fields are thought to have coherent magnetic 
fields with comparable strength to fluctuations. The major stumbling block to explaining 
these coherent magnetic fields by dynamo action in fully ionized gas is its tendency of gener- 
ating too strong fluctuations, which unfortunately inhibit the growth of a coherent (mean) 
component by back reaction (Lorentz force) when the strength of a m ean magnetic field is 
far below equipartition value - the so-c alled alpha quenching problem (ICattaneo fc Hughes 
19961 : iGruzinov fc Diamond! ^99^ . Il996h . 



It is, however, largely unknown whether and/or how backreaction constrains alpha ef- 
fect in weakly ionized medium, such as the galaxies, ISM, molecular clouds, etc, with am- 
bipolar drift (slippage between magnetic fields and the bulk of fluid (neutrals)). This is 
partly because almost all previous works on the effect of ambipolar drift invoked strong 
coupling approximation (the drift between ions and neutrals is balanced by Lorentz force 
due to sufficiently frequent collisions between the two), which makes ambipolar drift mainly 
act as a nonlinear diffusion. Thus, the ambipolar drift has been prima rily advocated as 



a means of enhancin g the effective diffusion rate over Ohmic value (e.g.. iMestel &: Spitzer 
19561 ; IZweibell Il988l ). It is also attributed to the fact that the important dynamic effect 
of fluctuat i ons (turbulence and Lor entz force back reaction) has often been neglected (c.f., 
Boss! 120001 ; iFatuzzo &: Adamsl |2002| ). Interestingly, these two factors come in together as 
strong coupling approximation is likely to break down on small scales (i.e., for fluctua tions) 
where Alfven frequency is larger than the ion- neutral collision frequency (IKiml 119971 ). The 
purpose of this Letter is to present a first self-consistent mean-field theory of the dynamo 
in weakly ionized gas, by incorporating these important synergistic effects of turbulence and 
back-reaction, without invoking strong coupling approximation. We shall demonstrate that 
ambipolar drift reduces alpha quenching, even overcoming it in certain (but extreme) cases. 

Before proceeding to a mean field theory of the dynamo, some insight into the effect of 
ambipolar drift can be gained by considering the problem of diffusion of a mean magnetic 
field in the two dimensions (2D). In the case of fully ionized gas, it is well known that the 



turbulent diffusion in 2D is severely reduced by back reaction (ICattaneo &: Vainshteinlll991 



Gruzinov fc Diamond! 119941 ) . In weakly ionized gas, the turbulent diffusion is still reduced 
below its kinematic value while ambipolar drift can increase the critical strengt h of a mea n 
magnetic field (above which the diffusion is reduced) by a factor of ^Jv in T n (IKiml 119971 ). 
Here, z/j n and r n are the ion-neutral collision frequency and correlation time of neutrals. 
Thus, ambipolar drift offers the possibility of dissipating a mean magnetic field at turbulent 
rate for sufficiently large Vi n r n . In this Letter, we also provide a self-consistent mean field 
theory for the diffusion of a mean mag n etic fi eld in 2D in weakly ionized gas, which not only 



confirms the numerical results of Kim 



(119971 ) but also generalizes the perturbation analysis 
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t herein. 



The remainder of the Letter is organized as follows. Section 2 contains the derivation 
of the effective dissipation rate of a mean magnetic field with ambipolar drift in 2D. We 
provide a mean field dynamo theory in weakly ionized gas in §3, by deriving the analytic 
expressions for alpha and beta effects. Concluding remarks are provided in §4. 

2. TURBULENT DIFFUSION OF MEAN MAGNETIC FIELD IN 2D 

In weakly ionized medium, Ohm's law is valid provided that the velocity of ions is used 
in the calculation of current. So, to consistently treat this problem, one needs to solve 
the momentum equation for ions and for neutrals (without Lorentz Force), together with 
induction equation for magnetic fields. Since v in j ' p n = v ni j ' pi and pi/p n <C 1 for weakly 
ionized medium, the neutral-ion collision frequency can be neglected. Here, p n and pi are 
the density of neutrals and ions, respectively, and v in and v ni are ion-neutral and neutral-ion 
collision frequency, respectively. Thus, the momentum equation for the neutrals is entirely 
decoupled from that of the ions as well as from the induction equation. We thus assume 
that neutral velocity is turbulent with a prescribed statistics, and then solve the momentum 
equation for ions, which evolves self-consistently by frictional coupling to neutrals and by 
Lorentz force. Note that we are not invoking the strong coupling approximation. 

In 2D, we work with ion vorticity u (uz = V x v) and magnetic potential A (B = 
V x (Az)), which are governed by the following set of two equations (in dimensionless form): 



Here, v is the ion velocity, N is the vorticity of neutrals, 7 = ^ n r n is the frictional coupling 
between ions and neutrals, 77 is Ohmic diffusivity, and \ is ion viscosity. We shall assume 
unity magnetic Prandtl number (i.e., rj — x)- By decomposing fields into large and small 
scale parts and assuming that there is no large scale displacement of the medium, we let 
v = v + v = v, to = ujq + uj = u, N = N + N = N, B = B + b, and A = A + a. 
Here, subscript '0' denotes a mean component, averaged over the statistics of N. Using this 
decomposition, we can separate the system eq. ([T]) in large and small scale components to 
obtain 



(d t + v • V)cj 
(d t + v • V)A 



-(B- V)V 
r]V 2 A. 



2 A + X V 2 u + -f(N -u), 



(1) 



[dt - XV 2 + 7 ]u; 
[d t - V V 2 ]A 
[d t - ??V 2 ]a 



(B - V)V 2 a + 7 Af, 
(v • Va) = -V ■ G , 
v • VA . 



(2) 
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Here, angular brackets denote the average over the statistics of N; in the first equation, the 
term V 2 A was neglected due to the large-scale variation of A Q ; in the last equation, the 
term v ■ Va — (v • Va) was dropped (quasi linear approximation). G = (va) is the flux of 
magnetic potential, which determines the evolution (effective diffusion) of A . To obtain G 
in terms of mean quantities, we rewrite it as: 

G = (v J dtd t a) + ( J dtd t wa) = Gi + G 2 . (3) 

Here, Gi is a kinematic part while G 2 comes from the back reaction of the flow onto 
the magnetic potential. It i s easy to check Gi = — ^(v 2 )VA , by using r approximation 



(IGruzinov &: Diamondlll994j . Il996l ). namely by replacing the time derivative by -. The ex- 
pression for Gi is just the standard beta effect in 2D. It is interesting to express Gi in 
terms of N since the statistics of the latter can be prescribed. For simplicity, we assume 
the statistics of N to be stationary with a delta-function power spectrum around k = ko as 
follows: 2 

(JV(k ls t)N(k 2 , t)) = ^5(ki + k 2 )5{h - k ) . (4) 

By taking spatial Fourier transform of the first equation of eq. ([2]) without the Lorentz force 
term, and by using eq. (JH), we obtain 



r 7 k 2 Q 



('" 2 ) = [T^rl 2 ^- (5) 



Therefore, the effective diffusion coefficient in the kinematic limit is given by f3 = = 

2"[ ]^~ ] 2 ^p^- Note that (3q takes its maximum value when neutrals and ions are strongly 
coupled with jt ~^>1. This is a natural consequence of the assumption that ions obtain their 
kinetic energy through frictional coupling to neutrals. Thus, crudely put, (3q is reduced by a 
factor [^] 2 . If there were an independent energy source for ions, this would no longer be 
true. 

To compute G 2 , we incorporate Lorentz force in the first equation of eq. fl2]) and take 
the Fourier transform to obtain: 

£(k) = —!—[-e ij3 [ dk'(k - k') i io(k - k')^A;' 2 S(k / ) + jN(k)} , (6) 

from which G 2 follows as: 

G 2 = -i j dkdk'e^ +k ^^e 3lm (u(kya(k')) (7) 

= ~ on I > V 2 a)VA --^- / dkdk'e^ +k ^^6 3lm (N(k)~a(k')) . 
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To calculate the first part of the RHS, we assumed that the magnetic potential fluctuations 
were isotropic and homogeneous. Since the neutrals are unlikely to be correlated with the 
magnetic field, the second part of the preceding equation can be neglected, leading to 

Compared to (3 = t(v 2 — b 2 )/2 in the fully ionized gas, the contribution from the backreaction 
in eq. (JE]) involves a multiplicative factor 1/(1 + tj). It is because the response of ions and 
magnetic field are different due to frictional coupling of ions to neutrals (see eq. d2]))- T o 



express (aV 2 a) in terms of large scale quantities, we use Zeldovich theorem ( jZeldovichlll957l ). 
which can be derived from the conservation of (A 2 ) in 2D-ideal MHD (by multiplying the 
third equation of eq. (J2J) by A and taking average over large scales) as 

v (aV 2 a) = (av) ■ VA = -(3(VA ) 2 . (9) 

Thus, from eqs. (|HJ) and (J9j), we obtain d t A — (rj + 8)V 2 Aq with 



8 = — . (10) 

Note that (tj + 8) is the total effective diffusivity of A Q . In the case of weak coupling limit 
(T7 C 1), the previous equation reduces to the beta suppression in fully ionized gas for a 
given p . Note, however, that (3q itself is proportional to (r7) 2 /(l + T7) 2 ~ (t7) 2 - I n the 



oppo site strong coupling limit (t7 ^> 1), one recovers an expression similar to that of (IKim 
19971 ) as (3 ~ 8q/(1 + B 2 /r] r y). Thus, back-reaction becomes insignificant when large scale 



magnetic field is weak enough so as to satisfy the following condition: 

(VA ) 2 = B 2 << rn = -^V> , (11) 

where R m = T(v 2 )/r] is the magnetic Reynolds number. Thus, the critical strength of mean 
magnetic field for the suppression of j3 effect is 7(f 2 ) / R m , which is larger by a factor of 7 
than that in the fully ionized gas. Note that the turbulent diffusivity can reach its kinematic 
value (3q = t(v 2 )/2 as 777 —» 00 but ca n never be greater. That is, (3 is the upper limit on 



(3. These results agree with iKiml (119971 ). 



3. MEAN FIELD DYNAMO THEORY IN 3D 



We now provide a mean field dynamo theory in weakly ionized gas in 3D, by self- 
consistently computing the alpha and beta effects. As previously, we use quasi-linear theory 
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to obtain the following set of non-dimensionalized equations for fluctuations and mean field 
(denoted by a subscript '0'): 

[<9 4 + 7-xV 2 ]v = B - Vb + b- VB - Vp + 7 N, (12) 
ld t - V V 2 }b = Vx (vxB ), 
[d t - r?V 2 ]B = Vx(vxb) = VxE. 

Here, N is neutral velocity (not vorticity); E = (v x b) = aBo — /3V x Bo is the electromotive 
force, which contains a and (3 effects. To compute E, we again consider two parts - the 
kinematic part E and the part coming from the back reaction of the magnetic field onto 
the fluids E x : 



E = (v x J dtd t v) + ( / dtd t v x b) (13) 
= a B - /3 V x B + (J dtd t v x b) = E + Ei , 



wher e a = — |(v ■ V x v) and /3q = |(i> 2 ) are the kinematic values (c.F., iKrause Radler 



19801 ). These two coefficients can again be expressed in terms of N as: 



«o = -^f(N.VxN), (14) 

ft = ^[TX-] 2 (iV 2 ). (15) 
3 I + T7 

The computation of Ei can most easily be done in Fourier space because of the pressure 
term. Thus, we write the equation for the velocity in Fourier space and then plug it into the 
electromotive force expression to obtain: 

Ei a = - ^ e a/3y J dp J dq Tp XlM (k - p)5 0A (q)(6 M (k - p - q)6 7 (p)) , (16) 

r a/ 3 7 (k) = 5 a pky + 5 ai kf3 — 2k a kf 3 k 1 /k 2 . 

To compute E\ a , we assume that the statistics of small scale magnetic fields is homoge- 
neous and isotropic but not necessarily invariant under plane reflection, with the following 
correlation function: 

(6 Q (k)6 )S (k')) = 5(k + k')[^-(5 aP - k a k p /k 2 ) + ip^e^kj = 5(k + k')^(k) , (17) 

where M(k) is the magnetic energy spectrum tensor and F(k) is the magnetic helicity spec- 
trum tensor. By using eq. ( 1TT|) in ( ITBl) and by keeping terms up to k (stretching and diffusion 
term), we obtain: 

E la = e aPy B 0X (k) y'dp$^(p)[(2^^-^-^p A ) + (5 /3A -^)^]. (18) 
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Since all integrals with odd numbers of pi vanish, E\ a reduces to: 

Eia = eafrBoxQ*)^ J dp ^^(^7 -P^/P 2 )(^x ~ ~^r) ( 19 ) 

f , ^(p) {0 Pf3P\P» x x M 

y p 8^p4 wp ^ ~ ^ — W/. - wi • 

The first part (proportional to fc^), contributing to /3, vanishes when integrated over angles, 
while the second part gives the correction term to a due to back reaction. Thus, there is no 
change in j3, namely, the turbulent diffusion is not affected by the back reaction of small scale 
magnetic fields in 3D with ambipolar drift! This result sharply contrasts to the claim made 
in the literature, based on strong coupling app roximation, that am bipolar drift enhances the 



diffusion of a me an magnetic field in 3D (e.g ., ISubramanianl Il998l ) . Note that the result for 



fully ionized gas (jGruzinov fc Diamondlll996l ) is recovered simply by taking the limit T7 — > 0, 
but by keeping (3 constant. On the other hand, the surviving part of Ei, contributing to a, 
reads: 

Ei a = — B 0a . (20) 

1 + r7 3 



Therefore, a effect, including the back reaction of the magnetic field, follows from eqs. ( fl"3l) . 
dHD, and (HDD as 

r (b • V x b) 
a = a + -— X - '-. 21 

1 + r7 3 

Note that only the helical (resp. non-helical) part of the magnetic spectrum is involved in 
the alpha (resp. beta) effect since a (resp. (3) is a pseudo-scalar (resp. scalar). Compared to 
the fully ionized gas, the contribution from the current helicity to a contains the additional 
multiplicative factor of 1/(1 +tj). This is again because the response time of ions is different 
from that of magnetic fields due to frictional coupling to neutrals (see eq. f|T2l ). Thus, it 
is very likely that the cancellation between fluid and current helicity for Alfven waves (as 
happens in fully ionized case with 7 = 0) may be avoided for jt > 1, thereby reducing the 
suppression of a effect. This shall be shown below. To close the expression for a, we need to 
express the current helicity in terms of mean magnetic fields. To do so, we use the topological 
invariant of mean magnetic helicity (a - b) in 3D, from which an analog of Zeldovich theorem 
can be derived as 

r)(b ■ V x b) = -(v x b) ■ B = -aB\ + f3 B ■ V x B . (22) 

Finally, combining eqs. (12T]) and (|22|) . we obtain the non- linear a-effect expression for 3D- 
MHD with ambipolar drift: 

«o+3^)B .VxB 

a= j 1 r Bl • ( 23 ) 

1 ^ 3(1+t 7 ) V 
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The previous equation recovers a in the case of fully ionized gas as 7r — ► 0. In the strong 
coupling limit (7T ^> 1), a effect is suppressed when eq. (fTTj) is satisfied. Therefore, the 
critical strength of mean magnetic field for the suppression of a effect is ^(v 2 ) / R m , larger 
by a factor of 7, compared to the case of fully ionized gas. This has significant implications 
for a mean field dynamo in the galaxy, ISM, etc where the bulk of fluid consists of neutrals 
with 7 ^> 1 (see §4 for more discussion). 



4. CONCLUSION 

We have presented self-consistent mean field theory of the turbulent diffusion (in 2D) 
and the dynamo (in 3D) in weakly ionized gas, by incorporating turbulence and back reaction 
of fluctuating magnetic fields. The key results are that in 3D, the backreaction does not alter 
the beta effect while it suppresses the alpha effect when the strength of a mean magnetic 
field exceeds the critical value B 2 ~ 7(t> 2 )/_R m . This critical value is larger than that 
B 2 ~ (v 2 ) / R m in the case of the fully ionized gas for 7 = v in r n > 1. Alternatively put, 
the suppression factor for the alpha effect is reduced by a factor of 7, compared to the fully 
ionized gas. The upper bound on a is given by its kinematic value ao = — r(v- (V x v)/3. In 
2D, the turbulent diffusion (/3 effect) was shown to be suppressed by back reaction when a 
mean magnetic field reaches the same critical value B 2 ~ j(v 2 ) / R m , with the upper bound on 
turbulent diff usion given by the kinematic value (3q = r(v 2 )/2. These results are consistent 



with those in iKiml (119971 ). 



Therefore, in weakly ionized gas, the degree of alpha quenching (in 3D) and the sup- 
pression of turbulent diffusion (in 2D) crucially depends on 7 = v in T n in addition to R m , 
i.e., the property of medium such as ionization, turbulence, etc. As v in ~ 10~ 2 n n cm 3 /yr 



(e.g., see 



Kiml 119971 ). v in r n ~ 10 n n for r n ~ 10 yr. Here, n n is the number density of neu- 
trals in unit of cm -3 . Therefore, in the limit of a very low ionization, the alpha quenching 
(and beta quenching in 2D) can be significantly reduced. For instance, in the case of young 
galaxy with n n ~ 1cm" 3 , L ~ lOOpc, v ~ lOkm/s, and T ~ 10 4 K, rj = 10 7 (T/10 4 )- 3/2 cm 2 /s 
~ 10 7 cm 2 /s, and R m ~ vL/rj ~ 10 19 . Thus, u in r n /R m ~ 10 5 /R m ~ 10" 14 , with the critical 
strength of mean field B c ~ 10~ 7 x \J (v 2 ), which is too weak. However, for dark molecular 
clouds with n n ~ 10 7 cm -3 , L ~ lpc, v ~ lkm/s, and T ~ 10K, rj = 3 x 10 11 cm 2 /s and 
R m ~ vL/rj ~ 10 12 . Thus, v in T n /R m ~ 1, giving B c ~ a/ (v 2 }\ Therefore, in this extremely 
low ionized gas, a mean field dynamo may work efficiently without alpha quenching. 

These results essentially come from the fact that the turbulence in weakly ionized gas 
does not become Alfvenic as the motion of ions undergoes frictional damping due to the 
coupling to neutrals (or, due to ambipolar drift). This is quite similar to what happens in a 
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very viscous fluid with R e <C R m {Re is the Reynolds number), in which case the suppre ssion 
facto r for alpha quenching is also reduced because of viscous damping of ion velocity ([Kim 



19991 ). Since R e <C R m i n the galaxies with a low ionization, the combined effect of ambipolar 



drift and viscous damping of fluid may render the mean dynamo sufficiently efficient, without 
severe alpha quenching. This interesting problem will be investigated in a future paper. 
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